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Abstract 

This paper is an appendix to a previous paper: |quant-ph / 010112^ "Relaxation Method 



for Calculating Quantum Entanglement", by Robert Tucci. For certain mixtures of 
Bell basis states, namely the Werner States, we use the theoretical machinery of our 
previous paper to derive algebraic formulas for: the pure and mixed minimization 
entanglements (i.e., Epure and Emixed), their optimal decompositions and their en- 
tanglement operators. This complements and corroborates some results that were 
obtained numerically but not algebraically in our previous paper. Some of the al- 
gebraic formulas presented here are new. Others were first derived using a different 
method by Bennett et al in |quant-ph / 9604024 . 



1 Introduction 



This paper is an appendix to a previous paper|jTl| by the same author. We will assume 
that the reader has read our previous paper. Without having done so, he/she won't 
be able to understand this paper beyond its Introduction. 

Henceforth, we will use "min." as an abbreviation for the word "minimiza- 
tion". In our previous paper we defined two quantum entanglement measures, the 
pure min. entanglement (Epure) and the mixed min. entanglement {Emixed)- These 
measures apply to any bipartite density matrix (the subscripts refer to the type of 
minimization space used, not to whether the density matrix is pure or mixed.) We 
showed that Ep^re is equal to the entanglement of formation, a measure of entan- 
glement first defined by Bennett et al in Ref.[^. Emixed-, on the other hand, is a 
new animal. We gave a numerical method for calculating Ep^re and Emixed-, their 
optimal decompositions, and also their entanglement operators (operators whose ex- 
pectation value gives the entanglement). We gave numerical results obtained with 
Causa Comun, a computer program that implements the ideas of Ref.[0]. We did this 
for a special type of Bell mixture called a Werner State and for Horodecki States that 
exhibit bound entanglement. 

In Ref.@, Bennett et al derived an explicit algebraic formula for the entangle- 
ment of formation of any Bell mixture. In Ref. 0], Wootters went one step further and 
generalized the formula of Ref . to encompass all density matrices of two qubits. 

In this paper, we use the theoretical machinery of our previous paper to derive 
certain algebraic formulas for Werner States. Specifically, we give explicit algebraic 
formulas for Epure and Emixed-, their optimal decompositions and their entanglement 
operators. This complements and corroborates some results that were obtained nu- 
merically but not algebraically in our previous paper. Most of our formulas for Epure 
were first derived, using a different method, by Bennett et al in Ref. 0. Our formulas 
for Emixed are new. 

2 Notation 

We assume the reader is familiar with the notation of Ref.|^. In this section we will 
introduce some additional notation that is used throughout this paper. 

We will use the notation of Ref.|^ intact except for one small modification. 
Ref. Jll dealt with a Hilbert space Tixy = Ti-x® 'Hy. Its two parts were represented by 
the random variables x and y (Xerxes and Yolanda). Here we will rename the two 
parts a, b (Alice and Bob). This conforms more closely with the rest of the literature. 
Also, it looks better in cases such as the one considered in this paper where one also 
uses X, y, z for indices of Pauli matrices. In conclusion, throughout this paper, we will 
be dealing with Hat where Sa = Sb = Bool. 

For any Hilbert space Ti. and any G Ti, we will often represent the projec- 
tion operator \ip){ip\ by 7!'{ip). C{7i) will denote the set of linear operators acting on 
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H. 

Let Zj^k be the set of integers from j to including both j and k. Let 
Bool = {0, 1}. Let x*"" be the n-tuple with x repeated n times. For example, 
0#3 = 0,0,0. 

The Kronecker delta function 6{x, y) equals one ii x = y and zero otherwise. 
We will often abbreviate 6{x,y) by 6y, 6{x,y)6{p,q) by 6^^, etc. Also, we will use 6'^ 
as an abbreviation for 1 — 6^,. In other words, is an an indicator function which 
equals 1 whenever jJ. ^ v and zero when jj = u. For example, if /i, z/ G ^0,3, then 
the metric in Special Relativity can be written as g^u = (5°° — One must be 

careful not to confuse S^^ = (1 - - 5^) and = 1 - S^S^. 

We will often use the color summation convention^. By this we mean that 
the summation signs will not be shown; summation will instead be indicated by 
displaying summed indices in a different color than the unsummed ones. For example, 
F^vV" = Ff_iuV'^- This is a better notation than the Einstein implicit summation 
convention which it is meant to replace. In the Einstein convention, we are instructed 
to sum over repeated indices. This becomes clumsy and requires a warning to the 
reader whenever we wish to use repeated indices that are not summed over. 

As is common in Relativity texts, we will often use Greek letters to represent 
indices that range over Zq ^ and Latin letters to represent indices that range over Zi 3. 
Unlike Relativity texts, we will not distinguish between upper and lower indices. 

Define JF2 = diag{l, —1,1). For any 3-dimensional vector n = (ni, n2, ns)-^, 
J^2n = (rii, — 77-2, ^3)^. One can likewise define for j G Zi^^ to be an operator that 
"flips" the jth component of the vector it acts on. 

For any 3-dimensional vectors a and x, a x x = Ax, where 
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a2 
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— ai 


— a2 


ai 






We will often represent A by [a x 
Let 





io)= , ii)= r (2) 



Let cr° = 1. Let a be the 3 dimensional vector of Pauli matrices. The Pauli matrices 
are defined by 



""^=[1 ) ^ ''y = [^ J ' = 1, -1 J • 

As is well known, the Pauli matrices satisfy: 

aV" = 61 + tekrjCT^ , (4) 
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for k,r,j G Zi ^, where ekrj is the totally antisymmetric tensor with ei23 = 1. Unfor- 
tunately, there is no formula for a^a'^ with /i, G Zq^s, that matches the conciseness 
and standardization of Eq.(H). Here is one particular attempt. 



where 
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(5) 



(6) 



and 



III/ 
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1111 

1 1 i -i 
1 -z 1 z 
1 I -i 1 



(7) 



Note that the operation n ® u defined by Eq.(|]) specifies an Abelian group 
(the operation is commutative, associative, has an identity, and has an inverse for 
each of its elements). The Abelian group defined by © on Zq.s can be shown to be 
simply the product of two copies of the group of two elements. 

Instead of defining f^^ by the table Eq.(|^), one can define it by the rather 
clumsy expression: 



f _ I \00 I rOO , rOO • 

f^y has a few useful properties. For example, it is Hermitian and it satisfies: 

For any x G [0, 1], the binary entropy function h{x) is defined by 
h{x) = — xlog2(x) — (1 — x) log2(l — x) . 



(8) 



(9) 



(10) 



Occasionally, we will also need to use h{x) with the base 2 logs replaced by base e 
ones. So define 



he{x) = (ln2)/i(x) = —xln{x) — {1 — x) ln(l — x) 



:iii 
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3 Bell Basis 



In this section we will discuss various properties of the Bell Basis. 

One can define operators that act only on the Tia (ditto, Tib) part of Tiab. Let 
A'^ (ditto, B^^) for /i G Zq^s represent the Pauli matrices that act on space Tia (ditto, 
TYb). Another natural notation for these operators is and a^. 

The following four states are usually called the "Bell basis" (with the magic 
phases) of T-igh- 

|i?(0)) = |=+) = ^(|00) + |ll)), (12) 
\B{1)) = zB'\ =+) = -^(|01) + |10)) = z| , (13) 
\B{2)) = zB'\ =+) = ^(|01) - |10)) = -I , (14) 

\B{3)) = zB'\ =+) = -^(|00) - 111)) = z\ =-) . (15) 

(By taking matrix products and linear combinations with real coefficients, of the 
operators 1, ia^, icTy and ia^, one generates what is called the Quaternion Algebra, 
invented by Hamilton.) 

The Bell basis states are an orthonormal basis of Tigh so they satisfy 

(i?(/.)|i?(^.)) = 5^ (16) 

and 

\B{f^)){B{^^)\ = l. (17) 

The Bell basis states place listeners a and b on equal footing: measurement 
of A'^ is the same as measurement of B^ up to a sign. Indeed, the action of A'^ on 
\B{0)) is 

A'^l =+) = {-lyiB'"] =+) . (18) 

The action of on |i?(z^)) for u ^ may have an additional —1 factor due to the 
fact that the Pauli matrices anticommute. For example, 

A^\B{3)) = AW\=+)=iB^A^\=+) = -'iB^B^\=+) = 

= BW\ =+) = B'^\B{3)) . (19) 

Thus, we see that, in general, the action of A^ on |-B(z^)) is 

Af'lBiu)) = {-lY'^i-lf>'Bf'\B{iy)) . (20) 
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Suppose fib (ditto, fig) is an operator acting only on Ht (ditto, TCg.)- We will 
call such operators local operators. It is easy to show that 

(i?(/i)|a|i?(/i)) = (=+ IQbl =+) = ^tr(fifc) , (21) 

and 

{B{ii)\na\B{fi)) = (=+ =+) = hiing) . (22) 

Note the right hand sides are independent of fi (although {B{fi)\flb\B{i')) generally 
does depend on /i and u.) Hence, all 4 Bell basis states contain the same amount of 
information about expected values of local operators. 

In future sections, we will need to find the matrix elements in the Bell basis 
of certain operators in CiTigh)- These operators can always be written as a linear 
combination x^i^A^B" . In this linear combination, will be acting on a Bell state 
so it can be replaced by plus or minus 5^. The product B^B'^ can itself be replaced by 
ffivB^^®'^. In this way, we can reduced the problem of calculating the matrix elements 
in the Bell basis of any operator in C{Ti.gb) to calculating the matrix elements in the 
Bell basis of B^. One has: 

{B{^,,)\B^\B{^,,)) = (-z)^^^/,,,eMJA..(0'-(=^|i?^®^^®^'l="') 

= (-0'-/A..(0'-5?e..e.2- (23) 
An equivalent way of writing the last equation is: 

{B{fi,)\B'\B{fi,)) = 6^^l , (24) 



and 



(fi(/il)|f ■5|5(/i2)) 



ioF 
-ix iix X 



(25) 



Another problem that we shall encounter in future sections is finding the par- 
tial trace with respect to either a or 6 of an operator X G C{T-Cg^). If 

X = x,,\B{^t)){B{u)\, (26) 

then one finds that 

traX = ^ {x^,^, + [xko - xok + XpgCpgk] } , (27) 

and 

trfeX = ^ {y^^,, + [yko - yok + Vpg^pgh] iA''} , (28) 
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where y^u = 1)''^ (— 1)*^'^. These equations generahze the well known result: 

trji?(/i))(i?(/i)| = i {n{\0),)+n{\l)t,)) , (29a) 

tn\Bi^^)){Bi^^)\ = i (TrdOy + 7r(|iy) . (29b) 

Thus, the Bell basis states do not distinguish between |0)a and (ditto, lO)^ and 
|l)fe) when the state of b (ditto, a) is unknown. An immediate consequence of the last 
equation is that the Bell states have maximum entanglement of formation. 

An element of C{T-Cab) can be expanded in various bases: one can expand 
it in terms of the operators A^B'^ for all /x, z/ G Zo,3 (call this the Pauli CiTiab)- 
basis), or the operators \a,b){a' ,b'\ for all a,b,a',b' G Bool (call this the Standard 
'C('^^ab)-basis), or the operators \B{fj,)){B{iy)\ for all jJi^u E Zq^^ (call this the Bell 
£(7Yab)-basis). In what follows, we will use mostly the Bell £(7iafe)-basis because it 
seems the most natural one for doing entanglement calculations. Thus, henceforth, 
whenever we represent C(Hab) operators by 4 by 4 matrices, the matrices should be 
understood as representations in the Bell C(Hab)-ha.sis. 



4 Entanglement of Pure State 

In this section we calculate the entanglement of any pure state of two qubits[0]. This 
is a good warm up exercise to prepare us for the following sections, where we address 
the harder problem of calculating entanglements of mixed states. 

Below, for any complex vector z, we will use \z\ = Vz ■ z* and z"^ = z- z. 

Any unit length \ip) G Hab can be expressed in the Bell basis as: 

\^) = {zo + iz-B)\=^) , (30) 
where = {zqI"^ + = 1. If 

p = |^)(^| = (^0 + tz- B)\ =+)(=+ |(^* - tz* ■ B) , (31) 



then 



where 



and 



tr„p = ^(zo + tz ■ B){z; - iz* ■ B) = no + n ■ B , (32) 



n = -{zqZ- zqz* + zx z*) . (34) 
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From Eq.(^) and Appendix 0, the eigenvalues of tr^p are simply no ± 

Hence, 



n\. 



Epure = Emixed = hijlQ + \n\) . (35) 

One can show using well known vector product identities that for any 4-tuple (^o, z) 
of complex numbers such that l^oP + = 1, one has 



Hence 



where 



zlz- zqz* + zxrY = l-\zf^ + zy . (36) 



^_.^_..(1±4HZ). (37) 



C=\zl + z''\. (38) 



C is called the concurrence^ of Ep^re is a monotonically nondecreasing 
function of C, and they both vanish at the same time, so C is also a good measure of 
entanglement. < C < 1. The pure state has C = 1 (maximum entanglement) 
iff its coefficients {zq^ z) are all real. 



5 Entanglement of Bell Mixture 

In this section we present the main calculation of this paper. For Werner states, 
we calculate Epure and Emixed, and their corresponding optimal decompositions and 
entanglement operators. Our calculation is split into 4 parts: (l)-R'"^, (2)/?"^ {?))Epure 

and Emixed (4)Aafe. 

We will call a Bell mixture any density matrix pab that can be expressed as 

Pat = T.^,\B{^)){B{^i)\ , (39) 
where J^^ii^fi = 1- We will call a Werner state any state that can be expressed as 

Pa, = mo\B{0)){B{0)\ + m,J2 \B{f^)){B{fi)\ , (40) 

where Dy G {1,2, 3} and mo + DyTUi = 1. This is a slight generalization from what 
is commonly called a Werner state. The term Werner state usually refers to the case 
where = 3 and the rank of pgb is 4. 
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5.1 K^fj Calculations 

We will begin by assuming that K^j^ can be expressed in a special "ansatz" form. 
Define 



h = dtag{l*'^\0 



(41) 



where Dy G Zi^^ was defined previously. Now assume can be expressed in the 
Bell representation as: 



-raT 



mo iqv 



(42) 



We assume that q and e are real, mo + D^mi = 1, and {{/"la G Zij^a} is a set of real 
3-dimensional vectors satisfying 



for all a G Zi^^a, 



a=l 



and 



(43a) 



(43b) 



a=l 



Here are some examples of sets of iT'^'s that satisfy Eqs.( 



= 2, iV, = 4, tTe 



-1 



J I, I 

V / \ 



/ i-ir 

{-If I \a,be Bool 

V 



(43c) 



(44a) 



(44b) 



= 3, iV, = 4, ^7- G <i I 1 I , I -1 



( 

1 

V -1 



(44c) 



Dy = 3,Na = 8, tT" G <( I (-1)'' \\a,b,ceBool 

-D' 



(44d) 
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Note that 



Pok = T.Kb = diagimo, (ml)#^^ O*^-^") 



(45) 



and 



Wn 



lab^^ =1/N^. (46) 

For Epurei it is clear that we want q = y/niomi and e = in Eq.(^2|). For Emixed, 
we intend to find those values of q and e that give the smallest possible conditional 
mutual information. 

Our ansatz i^"^ given by Eq.(^) depends on the following parameters: mo, 
mi, Z)j,, A'^, g, e and the v^. Out of these primitive parameters, one can construct 
the following auxiliary parameters whose use will significantly shorten our subsequent 
formulas. 



7] = D,- e{D, - 1) 



(47a) 



u = niQ + -qmi , 
mo — r]mi 



k 



Y=\q\jD, 



X = VWTy^ 



(47b) 
(47c) 
(47d) 
(47e) 



Next we will find the eigenvalues of K^f^. This can be done by using the 
following well known formula. Suppose M is a square matrix that can be partitioned 
into 4 blocks A, B, Ri, R2: 



M 



A i?2 

Ri B 



(48) 



where the submatrices A and B are square but Ri and R2 need not be. Then one 
can show that 



det(M) = det(v4) det(5 - i^i^iia) 



(49) 



One can use the last equation to find the eigenvalues of our ansatz K^i,. One finds 
(independent of a): 
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eigenvalue 


degeneracy 


K = 




/Na 


1 


A_ = 




/Na 


1 


Ao = emi/Na 


Du - 1 





3-D. 



(50) 



Note that since the eigenvalues of K^j^ must be non-negative, we must have 



e > and f - X > 0. Since 



-X' 



(51) 



it follows that |g| < ^mQimi^ < ^momi. One can also assume without loss of 
generality that q > since if g < 0, then one can replace q — g and -iT" —v" for 
all a. 

We also need to know InK^f,- To calculate In/T"^, it is not enough to find 
the eigenvalues of K^f^; we also need to find its eigenvectors. Our technique for 
finding them is inspired by Appendix 0, where we find the eigensystem of any 2 by 
2 Hermitian matrix. 

We begin by defining, for each a, three operators called E^, and P^^^: 



E, 



1 



tqv 



'aT 



p(0) 



-tqv^ 




/„ 



(52) 
(53) 
(54) 



Note that these operators satisfy the following multiplication table: 





Ea 




p{0) 

a 


Ea 


Ea 











Ea 





p(0) 

a 








P(0) 

a 



(55) 



Ea and Sq, can be used to define two new operators -P^^^-*: 

p(±) 



(56) 



Note that the Pjf") for a E = { — 1,0,1} satisfy the following multiplication 

table: 
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p{+) 

a 


p(-) 

a 


p(0) 

a 


p(+) 

a 


p{+) 
a 








p(-) 
a 










P(0) 

a 








p(0) 



(57) 



Thus, the P^'") are orthogonal projection operators. 

It is easy to show using the definitions of Pjf^ and Ao- for a E that 



E 

t6Z_i,i 



Thus, 



lni^^= Y: HK)Pt 

o-G^_l,l 



(58) 
(59) 



Technically, we should also add a term ln{0)diag{0'^^''~^^ , to the right hand 

side of the last equation to account for the 3 — zero eigenvalues of K^i,. However, 
we can safely ignore this infinite summand if we only use In /T"^ in expressions where it 
is multiplied times Pab. The infinite summand is annihilated when In K^f^ is multiplied 
times Pab- 

5.2 i?"^ Calculations 

To find R^h, we need to calculate the partial traces of K^f^. One gets 

1 n 



tlaK, 



ab 



V2 



■ B 



and 



where 



J^a — ^'-bJ^ab 



1^ + (^2n°) ■ A 



^ No, V2 



(60) 



(61) 



Therefore, 



n 



gir 



p: 



ab 



Wo, 



V2 



- + (J-^n") ■ A 



We also need to know InP"^. Using Appendix 0, one finds 



(62) 



(63) 



InP^ = -lniV« + lnQ + n"-p)+ln(^i + (J-2n")-i) 

= -lniV^+ In f ^ + (-1)^>^) n , (64) 



^eBooi 
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where 



(65) 



for ^ G Bool. 

At this point we have calculated In i?"^, but we have not yet expressed it in the 
desired form, as a matrix in the Bell representation. To do this, we need to find the 
matrix elements in the Bell basis of the projectors vr(|^^)b) and T^{\^r)b) for ^ G Bool. 
These matrix elements can be found using the techniques discussed in Section ^. One 
finds: 



1 r 

2 

1 



1 



-ifB 



+ 



r 

ir^ 
-if i(fx 



and 



vr(le-)a) 



1 + (-1)^1- f 







Putting all this together, we get 



(66) 



(67) 



In R 



,1 



- -IniV^ + ln ( (^ + r)(- - F)j +ln f-^ 



m""^ 
-m" 



(68) 



5.3 £pme and Emixed Calculations 

Recall from Ref. [ll] that the following Lagrangian C must be minimized to obtain 

both Epy^rf. and -Emixed- 

C = Ik-Ir, (69) 

where 



/i. = Etra6(i^a6lni^."J 



(70) 



and 



lR = Y.iT^iK^^^R"J ■ 

a 

Using the results of previous sections, one finds 



(71) 
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lK = -\nN^+Y.X'^ ln(A;) + {D, - 1)X', In(A'o) , (72) 

where A'^ = NaXa for a e Z^i^i, and the X„ are just the eigenvalues of K^f^ that we 
found earher. One also finds that 

ln = -\nN^-2K{^ + Y) . (73) 
Putting all this together, we get 

f E.=±(| + crX)ln(| + aX) 
C = l +(D^ - l)emiln(emi) . (74) 

[ +2K{l + Y) 

Next we will use Eq.(|T^ to calculate entanglement E = mm(£)/(2 ln(2)) for pure 
and mixed minimizations, 
(case l)Pure Min. 

In the case of pure minimization, one has q = ^/rnoffTi and e = 0. Thus, 
the auxiliary parameters defined by Eqs. (^7[) reduce to: r] = D^, u = 1, k = {uiq — 
D^mi)/2, Y = ^mQ^D^mi) = ^uiq^I — mo), and X = 1/2. Thus, from Eg . ([TiD , we 
get 

Ep^re = h(l + Y) . (75) 



.2 

If we define the concurrence C for this case to be: 

C=|2mo-l|, (76) 
then Eq.([75|) can be rewritten as in Ref.[0: 



E,^re = h (^^4^) • ^^^^ 

(case 2) Mixed Min. 

For mixed minimization, the constraints q = ^/rn^ymi and e = are no longer 
required in order to make i^"^ separable. We can choose q and e so as to minimize C 
given by Eq. ([7^) . Treating £ as a function of q and e and setting its partials to zero, 
we get the following two constraints 

de \ or mi{D^ - 1) = ^ V ) 



and 



f = ^ ^ln(^-^)=^ln(^^) . (78b) 

















if) 
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Figure 1: Plot of C versus Y at fixed k when D^, = \. See Eq. (|80D . 



In general, 

^ mixed TTi o ' C^"^) 

/m/ 

with £ given by Eq.(^), subject to the constraints Eqs.(|78D. 

Eqs.([78|) with ^ \ constitute a system of two nonlinear equations which 
one would like to solve for the two unknowns e and q. As far as we know, this "e — q 
System" cannot be solved exactly in closed form- its roots can only be found using 
numerical techniques. See Appendix ^ if interested in an approximate analytical 
solution of the system. 

The case of mixed min simplifies considerably when = 1. In that case /T"^ 
and L are independent of e, so we need only minimize over q. L reduces to 

c = -K + v^FTT^) + 2K + ^) • (80) 

Fig.|l| is a plot of C versus Y at fixed k, according to Eq. (|80|) . The largest possible 
Y value corresponds to pure minimization. £ does not achieve its minimum at that 
endpoint, but rather at a smaller value of Y , which we call Yq in Fig.|l]. l^o can be 
determined by solving Eq.( [78b|) for F as a function of k. We conclude that for D„ = 1, 



Emixed = -(^)^(^ + \/k' + Yi) + + Fo) • (81) 
5.4 Calculations 

Next, let us calculate for the cases of mixed and pure minimizations, 
(case l)Pure Min. 

The entanglement operator for pure min satisfies: 

ACTV^a) = (lnK:,-ln<,)|V;.), (82) 
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where 



/rrio 



3) 



—ty/miv I 

Note that {ipali^a) = 1- Using the results of previous sections, one can express the 



right hand side of Eg . (1821) as follows. Define auxiliary quantities /, M and a by 

/ = 



mo 



M 



In 



diag 



J 



and 



Then 



Since 



we get 



1 1 
In ( (- + mo/)(- - mo/) 



ln<,|7/.„) = (-lniV^-M-a)|^„) 



^5) 

^6) 
^7) 



A^- = M + a. (89) 

Using this value for A^J^^*^ and the value for Epure that we obtained in Section p73| , 
and also using the constraints g = ^ym^rn^, e = 0, one can check that 



[2 In 2) Epure = tr(pabA 



pure\ 
ab ) 



(90) 



(case 2) Mixed Min. 

The entanglement operator for mixed min satisfies: 



A mixed 



lnJC,-lni?:,. 



(91) 

Using the results of previous sections, one can express the right hand side of Eq. (pl|) 
as follows. Define the auxiliary quantities M and a by 



M = diag 



E [I + ^) In + -X) , (ln(emO)#^% 0#3-^^ 



(92) 
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a = -ln((l + F)(l-F)) . (93) 

Then 

^nnxed = M + a. (94) 

Using this value for A^"'*^*^, and the value for Emixed that we obtained in Section [5.3| , 
and also using the constraints Eqs.(|78D, one can check that 



2 In 2)E^,^,, = tr(p^A™^'='^) . (95) 



6 Implications 

If f2 is an orthogonal matrix and we replace our vectors t?" by fit?", then K^^, changes 
but the value of the entanglement doesn't. Thus, there is a continuum of possible K^^^s 
that minimizes C{K,K). It is convenient to define a pure-min- entanglement orbit 
(ditto, mixed-min- entanglement orbit) as a set of all K^h G JCpure (ditto, K^^^ G JCmixed) 
which are stationary points of C{K, K), and which give the same value for £{K, K). 
Only one pure-min (ditto, mixed-min) orbit is a global minimum of C{K,K). Any 
K^fj (or any orbit) whose Aab is independent of a, will be said to be a-insensitive (or 
just insensitive for short). 

In previous sections, we found a pure-min insensitive A'"^ whose concurrence is 
|2mo — 1|, regardless of whether niQ > | or not. On the other hand, in Ref.[^, Bennett 
et al found a different K^f^ whose concurrence is zero when mo < ^. (In fact, Ref.0 
shows that any Bell Mixture, not just the Werner states that we are considering here, 
must have zero entanglement of formation when the largest weight mo < |). Thus, 
for the Werner states that we are considering here, there exist at least two pure-min 
insensitive orbits when mo < |. 

In Ref. [H, we claimed that there is only one orbit, the global minimum of 
C{K,K), that is insensitive. This paper has given a counterexample to that claim. 
It appears from this paper that the insensitivity condition is necessary for the global 
minimum orbit, but it is not always sufficient. This leads one to wonder why and 
when the sufficiency part of the proof in Ref.|jl| breaks down. The breakdown may 
be due to the fact that the proof treats in a cavalier manner the infinities produced 
by taking the log of zero eigenvalues. 

If we generalize our ansatz for K^f^ so that the -u^'s can be complex, can we 
find any more insensitive orbits? We try to answer this question in Appendices 
and 0. We find that going from real to complex iT°'s yields new insensitive orbits in 
the pure min but not in the mixed min cases. We also find that for the iP-complex 
ansatz, there exist a countable number of pure-min insensitive orbits, and each of 
these corresponds to a stationary point of the pre-concurrence. 
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A Appendix: Eigensystem of 2 Dimensional 
Hermitian matrix 



Consider any 2 by 2 Hermitian matrix n. One can always express ii a.s n = riQ + n- a, 
where uq and n are real. The eigensystem of n follows immediately from the following 
easily proven identity: 

no + rI-a=K + v^)f^^^l + K-v^)f^^^l ' (96) 



where h = nj Define 



Then 



and 



2 



p. = . (97) 



P+P_ = P_P+ = , (9^ 



{P±f = P± . (99) 

Thus, P+ and P_ are the projectors onto the two eigenspaces of n with respective 
eigenvalues Uq + and uq — V^. 

An alternative, more tedious way of finding the eigensystem of h is to rotate 
the equations (T^IO) = |0) and = — Define a rotation vector 6 by: 

9 = 9^:^, ^ = arccos^, (100) 
\z X ril \n\ 

The spin up and down states along the n direction can be obtained in terms of those 
along the z by: 

|0^) =exp(-z^)|0) , (101) 



2 



and 



|l^)=exp(-z^)|l). (102) 
One can show that the projectors P± defined by Eq.(^) satisfy: 

P+ = |0,,)(0^| , (103) 

and 

P-=|ln)(ln|. (104) 
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B Appendix: Approximate Solution of e—q System 



In this appendix, we will assume > 1. Eqs.(|78D constitute a system of two 
nonlinear equations which one would like to solve for the two unknowns e and q. 
As far as we know, this "e — q System" cannot be solved exactly in closed form- its 
roots can only be found using numerical techniques. Nevertheless, as we will show 
next, it is possible to get approximate analytical expressions for its roots. 

It is easy to show that the e — q System is solved exactly by (e, g) = (1,0). 
This root, however, does not yield the global minimum mixed- min orbit, so it is of 
little interest to us. It can be rejected if we restrict our attention to roots for which 
the regulator q is nonzero. 

From the numerical results obtained with Causa Comun and discussed in 
Ref. we expect that the e — q System has a second root which is very close to 
the pure min case, for which (e, g) = (0, y/mQmi). The rest of this appendix will 
be devoted to finding an approximate value for this second root. Taylor expansion 
at (e, q) = (0, y/monii) is not possible since both equations of the e — q System are 
non- analytic at that point. Another type of approximation is called for. 

We begin by rewriting the e — q System in the following equivalent form: 

/U \22X/U \ 2^2X 



xV 




(105b) 



When e ^ and q ~ ^Jra^tfil, one has that u^u^ = l^k^k^ = uiq — i, 

1 

2' 

Eqs. (|105| ) by the following two equations: 



y ~ yo = -i/mo(l — mo), and X X^ = ^. Our approximation consists of replacing 




u \' /I Y \ 



In effect, what we are doing is setting e = and q = ^Jln^ml everywhere in Eqs.( |105| ) 
except where a very small number (namely, | — X and \ — Y^ appears raised to a 
power. Since ^ — F° = iff mo = |, this approximation is expected to work best in 
the vicinity of mg = \. Define the quantity p by: 

p = m^mrr]-q^D., = {-f -X\ (107) 
Eqs.( |10ti| ) can be expressed in terms of p as follows: 



19 



(108a) 



.0,1 



1 

2 



qJD^ . 



(108b) 



Eqs.([T08|) imply 



momi?7 = momi [D^ — e^D^ — 1)] = momi 



A, 



mi 



(109a) 



and 



(109b) 



Now the left hand sides of Eqs.( |109| ) are just the two terms whose difference defines 
p. Therefore, 



0\2 



p = (Y^) 



P 



mo 



mi Dy 



2 ^2 



(110) 



Motivated by the last equation, we define a function / of p by: 



0\2 



/(p) = -p+(n 



p 



1 - mo 



{D, - 1] 



2 ^2 . 



(Ill) 



Eqs. (|108|) can rewritten as follows, so that they express e and q in terms of p: 



'112a) 



mo 



2-^ 



4 + 



(112b) 



Clearly, given any proot for which f{proot) = , one can use Eqs.( |112|) to calculate a 
point (e, q)root that approximately satisfies the original e — q System. If (e, q)root yields 
values for u, k, X and Y that are close to u^, /c°, X'^ and Y'^, respectively, then our 
original assumptions are vindicated and we say the approximation is self consistent. 

Let mo = I + ^fn- We wrote a simple Excel spreadsheet in which we plotted 
/(p) as a function of p, with inputs mo and Dy. Fig.^ shows schematically what we 
found. For Dy = 2, /(p) has a zero regardless of the sign of 5m. For Dy = 3, /(p) 
has a zero for 6m > but not for 5m < 0. 

We also calculated analytical approximations for /, and for its first and second 
derivatives, in the limit of small p and small \6m\: 



/(O) = {Y 



0\2 



1 



< 



(113) 
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D = 3 


8m>0 


f(p)i 




f(p). 


















5 m<0 


f(p)i 






f(p). 








p 



Figure 2: Schematic plots of /(p). Four cases depicted depending on whether D 
equals 2 or 3 and whether (5m is greater or smaller than zero. 



dp 



1 

7p 



5m 



-{Dy-l) 



- \5m\ 



if 5m < 



(114) 



3 

P2 



^(D„-l)(<5m-l)p^™-l 



^ if (5m > 

if 5m <0 



8pf + l*™l 



(115) 



C Appendix: Stationary Points of 
the Pre-concurrence 

Define the pre-concurrence amplitude by 



j=0 



(116) 



where nij > for all j G ^o,n; J2]=o^j = 1) the are real numbers, and we fix 
6*0 = 0. Next define the pre-concurrence Cg by 



Co 



;ii7) 



The global minimum of Cg (at n = 3) over all phases ^ is called the concurrence. This 
global minimum arose in our calculation of the pure min entanglement of a two qubit 
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system. One wonders if Ce has other stationary points, and whether they play a role 
in our theory. In this appendix, we will find the stationary points of Cq. In Appendix 
0, we will show that such stationary points are indeed very relevant to our theory. 




Figure 3: 751 is a "fractured unit vector" . 




Figure 4: min \'Jq\ equals zero if mo < | and mo — Sj=i if "^0 > |- 

From Fig.|^, we see that 7 is a "fractured unit vector" in the complex plane: it 
equals the vector sum of segments whose lengths ruj add to one. From this geometrical 
picture, the smallest possible length for 7^ is clear. As shown in Fig.^ if mo < |, 
then min|79| = 0. If mo > |, then min |76i| = mo — I]j=i"^i- Are there any other 
stationary points of Cg7 Yes. 

Any stationary point 6 of Cg must satisfy: 

n 

5Cl = Si^eYo) = ^ E SO^m^ie^'^Ye - e-^'^le) = . (118) 

Hence, 

m,|7,|sin(^,-Z(7e)) = 0, (119) 

for all j G Zi n- The last equation is satisfied iff either 7^ = 0, or, for all j G „ 
such that rrij 7^ 0, 6j = Z(76)) + irbj + 27mj, where bj G Bool and the rij are integers. 
Hence, the set of stationary values which Cg can assume is given by 
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^(-l)^^mjl : bj e Bool, bo = o\ , (120) 
i=o J 

together with Cq = 0, which, however, is only possible when mo < |. 
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Figure 5: Values of CeiOi, 62) shown in color. Cq is periodic in 61 and 62- Its minimum 
value is |mo — 2mi| if mo > | and zero otherwise. = (possible only if mg < \) 
is indicated by an asterisk. 

Fig.^ shows a plot of 

Cg{e^,e2) = |mo + mi(e^^^ +e^^^)| ; (121) 
that is, Cg when n = 2 and mi = m2. 

D Appendix: Complex i/"'s 

Recall that an orbit is a set of all K"f^ which give a common stationary value for 
C{K,K). And if Aat for an orbit is independent of a, we say that the orbit is 
insensitive. As was pointed out in Section |, the results of this paper imply that 
there can be more than one pure-min insensitive orbit. If we generalize our ansatz for 
K^i, so that the ?T"'s can be complex, can we find any more insensitive orbits? We know 
that the global minimum of the pre-concurrence corresponds to an insensitive orbit. 
Do the other stationary points of the pre-concurrence also correspond to insensitive 
orbits? In this appendix we will show that going from real to complex -u^'s yields new 
insensitive orbits in the pure min but not in the mixed min cases. We will also show 
that for the u^-complex ansatz, there exist a countable number of pure-min insensitive 
orbits, and each of these corresponds to a stationary point of the pre-concurrence. 
Let f/" be defined by 

= diag{e"^",e''''^,e'^^) , (122) 
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where the are real. Let 



1 

U"" 



(123) 



One can define complex vectors in terms the real ones -u" by: 

= U'^v" . (124) 



Note that by virtue of Eqs.(|D, C"'^C" = D^, EaC" = 0, and EaC"C""^ = NJ^. Our 
new ansatz for K^j^ is defined in terms of the old one by: 

= U^K'^jU^r = (K^),^^^^^ . (125) 

Henceforth, as we have just done, when we need to distinguish between a new quantity 
and the corresponding old one, we will indicate the new one by a tilde. (By old 
we mean with the iP real, and by new, with the complex.) Sometimes, if it 
is clear from the context that we are speaking of new quantities, we will drop the 
tildes. Alternatively, sometimes we will introduce a new symbol for a new quantity 
to distinguish it from the corresponding old one, as we did by introducing the symbol 
to represent the new -u". 

It is convenient to define zq and by: 

zo = y/m^ , (126) 
F = -i^C ■ (127) 



Note that {zoY + ^ ■ z"* = 1. 

One can define a concurrence amplitude 7 by 

3 

= mo-mi^e^'^^"(wj')2 . (128) 

The concurrence C is defined as the absolute value of 7: 

C=\l\. (129) 

Note that in general, 7 and C depend on a. However, as will soon become apparent, 
they must be independent of a for any insensitive orbit. 

Next we will calculate the new i^^^, and Ag^. The calculation is very 
similar to that of the old K"^, and A^b presented in previous sections, 
(a) K^i^ Calculations 
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Since K^j^ and K^t, are simply related by a unitary transformation, it is clear 
from previous results for K^f, (see Section that 



where 



lni^^=-ln7V^+ J2 HK)PL''\ 

cre^-1,1 



Written more explicitly, this becomes 



(130) 



:i3ii 



-lnAr„ 

+ E.=±ln(f + aX)| 

/ 
+ ln(emi) „ ^ ( 



1 1 ak 




2 ^ %X 


2X 






2X 













(132) 



(b) R^j, Calculations 

that: 



We proceed as we did in Section p.2| where we calculated i?"^. Now we find 

(133) 



1 /I 



^ No, V2 



= tiaK, = — - + (n" + An") ■ 5 , 



and 



where 



and 



kl = tnk^ = ^{l + [-^2(n" - An")] ■ I 



y/moTTli 



An" = -(i« X ^*){l-e) . 



(134) 



(135) 



(136) 



Note that the An" enters with opposite signs in Eqs. ( |133]) and ( |134| ). Proceeding as 



we did in Section 5.2, we finally find 



lni?:,= <! 



-IniV^ 

+ ln((l + F)(i-F) 



+ ln I 



h+y\ 1 



Y / Y 



:taT 





—in'^ iAn°' x ■ 



(137) 
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where 



F = |f?" + An"| . (138) 

Note that in general, Y depends on a, However, as will soon become apparent, 
Y must be independent of a for any insensitive orbit. Note also that • An" = so 
Y"^ = |n°p + |An°p. After some algebra, one can show that: 

|.lT=f^)^(l-7.), (139) 

and 



I An' 




2 



4 




4 



{(l-7r)[7r.-(2mo-l)]-7-} , (140) 



where jr and 7^ are the real and imaginary parts of 7. We see that when 7 = 1 for all 
a, n" = An" = for all a. In the other extreme, when 7 = 0, one gets |n°p = 

and |An°|2 = 11^(1 -2mo). If besides 7 = 0, one assumes the pure min constraints 
q = ^momi, e = 0, then |n"p = mo/2, |An"p = (1 — 2mo)/4, so y = |n" + An"| = 1. 
Recall that for the real-'u" ansatz, one has 7 = |2mo — 1| and Y = y^mo(l — nig), and 
therefore, 7 = implies mg = Y = ^. We see that in the complex- w" ansatz, it is 
possible to have 7 = 0, F = |, and mo 7^ |. 
(c) Aab Calculations 

For the mixed min case, A^"'^'^ is simply 

Ar^'^ = ln^^-ln^^, (141) 

where Ini^"^ and lni?J^^ have just been calculated. 

For the pure min case, (e, g) = (0, y/monii) and \ijja) is defined by Eg. (p3|) , 
which can be expressed in terms of zq and ^ as 

\^a)=( Z) ■ (142) 



z 



One finds 



Krii^o.) = (lni^^-ln^^)|^„) 
-ln((i + F)(i- 




\ }m- (143) 

l + ^(f/")2* I 
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As pointed out earlier in this appendix, if g = y/momi, e = 0, and 7 = 0, then Y = ^■ 
In this hmit, Eg. ( |143|) gives A^"''^ = 0, as expected. 

A simple Lemma: Consider an insensitive orbit for which ra" and An" are 
independent of a. If this premise is satisfied, then q = and e = 1. The premise is 
true in particular when 7 = 1 for all a. Proof: By Eq. (|137|) and the premise, Ini?^^ 
is independent of a. Since Ini?^^ and Aab are independent of a, In/C"^ (and therefore 
K"f^) must be independent of a too. From the ansatz form for i^"^, it is clear that if 
K^fj is independent of a, then q = and e = 1. If 7 = 1 for all a, then by Eqs.( |139| ) 
and (11401) , n"' = An"' = 0. Therefore, n"' and An" are independent of a. QED 

We end this section by finding all insensitive orbits that are possible according 
to the just derived formulas for InK"^, Ini?"^ and Aab- 
Pure Min Case: 

Our arguments are based on Eq. (|143|) , which we restate here for convenience: 



X pure 



{Ink, 



a 
ab 



lni?:,)|^.) 





m . (144) 



Zi,3, where hj G Bool is independent of a. This gives 7 



"^nij, where raj 



In the examples given by Eqs. 



We will henceforth call this equation Eq.A. Note the following. 

For the right hand side of Eq.A to be independent of a (for an open set of rriQ 
values), Y and 7 must be independent of a. 

If 7 = 0, then we must have rriQ < ^, as explained in Appendix y. According to 
that Appendix, 7 = is one of the possible stationary values of the pre-concurrence. 

If 7 7^ 0, then from Eq.A, ([/")^ must be real and independent of a. Thus, 

g*20- = for j e 
mo-E?=i(-i)'^ 

is independent of a for all j. 

If bj = 1 for all j, then 7 = 1. By the Lemma proven previously, this implies 
that q = and e = 1. We can exclude this case because q = and e = 1 would not 
give (except for some special values of mo) a i^"^ of the form Wa\'ilJ°'){ip°'\- 

Comparing the above results with Appendix y, we see that there are a count- 
able number of pure-min insensitive orbits, and each of these corresponds to a sta- 
tionary point of the pre-concurrence. 
Mixed Min Case: 

Our arguments are based on Eqs. (|132|) , Eqs. (|137|) and Eqs. (|141|) , which we 
restate here for convenience: 
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InK. 



ab 



2X 



+ ln(em 








In 
In 



A mixed 



Y / y 







at 



2JC 



m 

iAn" X ■ 



:i45a) 



(145b) 



(145c) 



Henceforth, we will call these 3 equations A, B and C, respectively. 

Subtracting entries (1,0), (2,0) and (3,0) of the right hand sides of Eqs.A and 
B, we get: 



In 



~2X 



-In 




(146) 



This vector must be independent of a. Since n° = |(C" + C"*)? ^6 must have AiC" + 
= V, where Ai, A2 are real numbers, and Ai, A2, w are independent of a. But then 
Ea(AiC° + A2C*) = = so A1C + A2C* = 0. In other words, (Ai + Ase-'^*^^")?;^" = 
for all a and j. Assume Ai and A2 are non-zero for an open set of mo values. Also 
assume 
in Eqs.( 



assume that 7^ for all a and j G ^i,d„ (this is true for the examples given 



. Then e~*^'^j must be real and independent of a and j G Zi^d^. If 
1, then 7 = 1. But by the Lemma proven previously, this implies that 
g = and e = 1. We can exclude this case if we require the regulator q to be non- 
zero. Thus, there appears to be just one mixed-min insensitive orbit, the one found 
in Section ^ for the real-t?" ansatz. 
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